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Abstract 

Collective coherent scattering of laser light induces strong light forces between polarizable point 
particles. These dipole forces are strongly enhanced in magnitude and distance within the field of 
an optical waveguide so that at low temperature the particles self-order in strongly bound regular 
patterns. The stationary configurations typically exhibit super-radiant scattering with strong par¬ 
ticle and light confinement. Here we study collective excitations of such self-consistent crystalline 
particle-light structures as function of particle number and pump strength. Multiple scattering and 
absorption modify the collective particle-field eigenfrequencies and create eigenmodes of surpris¬ 
ingly complex nature. For larger arrays this often leads to dynamical instabilities and disintegration 
of the structures even if additional damping is present. 
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I. INTRODUCTION 


Laser light scattering off free particles is accompanied by momentum transfer and thus 
forces. If the incident light helds are far detuned from any optical resonance, energy absorp¬ 
tion and spontaneous re-emission play only a minor role and a coherent redistribution of 
photon momenta determines the force. As such coherent scattering from different particles 
largely preserves the laser phase, these helds interfere and thus the resulting force on each 
particle depends on its distance to nearby particles and acquires a collective nature, which 
depends on all other particle positions [1,2]. While in a 3D geometry the interaction strongly 
decays with distance and is averaged out by phase randomization from particle positions 
and motion, in regular particle arrays or in a spatially conhned geometry, this interaction 
can become very strong [3, 4]. Particularly strong effects appear if the light helds are op¬ 
tically conhned in resonators [5], guided by optical structures as mirrors [6] or conhned in 
waveguides [7-10]. Besides cold dilute gases, closely related ehects have been studied using 
various kinds of nano-particles in solution [11-13]. 

Following hrst predictions by Chang and coworkers for near resonant weak pumping of 
atomic dipoles [7,10], we have recently shown that cold particles illuminated far of resonance, 
who can freely move along or within an optical nano-hber held, tend to form regular but 
complex optical structures through optical long range coupling. They collectively scatter 
light into the hber mode and self-trap in the optical potential generated by interference of 
the scattered light and the pump held [9]. Surprisingly we see that not only the light conhnes 
the particle motion but the particles also conhne the light forming Bragg like atom mirrors 
at their outer edges. Numerical simulations predict that in certain cases even more complex 
structures as cavity arrays are self-formed by the particles [8]. This potentially generates a 
self-ordered cavity QED system with strong particle photon coupling via so called atomic 
mirrors [14]. While such dipole energy minimizing conhgurations were also studied for 
conventional ID optical lattices, they could be shown to be generally unstable [15]. 

In this work we theoretically study the dynamics and stability of such self-organized cavity 
QED systems. Already in our previous work we could identify conhgurations where all forces 
on the particles vanish and one gets a restoring force towards these equilibrium positions 
for small shifts for each particle. Nevertheless, displacing one particle immediately changes 
the helds at all other particle positions and creates perturbation forces on them. Here we 
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cannot expect any energy conservation for the particle motion as we are dealing with an open 
system where the pump laser forms a non depleting energy and momentum reservoir [16, 17]. 
Hence perturbing one particle induces nonconservative collective motion of all particles. 
Eventually these oscillations exhibit exponential growth and lead to a disintegration of the 
whole structure. A central aim of this work is to calculate and study the frequencies and 
stability of the corresponding linearized eigenmodes of the coupled atom field systems and 
their application to induce tailored long range interactions in this system. Of course the 
system is inherently nonlinear so that the full dynamics can only be captured numerically. 

Our work is organized as follows. First we quickly overview the results of our recent 
paper [9]. Out of these, in order to investigate the stability of the particles, we linearize the 
forces acting on the particles around their equilibrium positions. Calculating the eigenvalues 
of the so obtained coupling matrix we hud a condition which enables and limits the formation 
of stable conhgurations. With this we first study configurations in the negligible coupling 
limit and then pass over to more realistic complex coupling parameters. 


II. SCATTERING MATRIX APPROACH TO LIGHT INDUCED MOTION IN ID 

Let us hrst briefly recall our previously developed multiple scattering approach to calcu¬ 
late the stationary helds and forces on a linear array of particles with transverse illumination 
as schematically depicted in Fig. 1 [9] . For a given spatial arrangement of point scatter- 
ers in the held of a ID waveguide the light held can be calculated by a scattering matrix 
approach [17-20]. The contribution of each particle is represented by a single 3x3 matrix 
parametrized by an ehective coupling constant ( proportional to its linear polarizability a 
and the held mode amplitude at its position. This determines the interaction between the 
particle and the hber mode as well as the scattering amplitude in and out of the hber rj. The 
transmission and rehection coefficients are related to the coupling constant via t = — 

and r = fC/(l ~ *0- symmetric scattering the coupling between the amplitudes left and 
right of a particle can then be expressed using the following beam splitter matrix [9] 
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Figure 1: A ID array of point particles scattering light in and out of an optical nano-structure can 
be modelled as a collection of beam splitters interacting with a plane wave. 

where Bj, Cj, rj are the incoming helds allowing to determine the outgoing fields Aj, Dj, cf. 
Fig. 1. These in turn are used as inputs to neighbouring scatterers as Cj-i = Aj exp{ik{xj — 
Xj-i)) and Dj_i = Bj exp{—ik{xj — Xj-i)). Multiplying all the matrices including the 
propagation phase shifts then allows to determine the held distribution corresponding to 
the momentary particle positions. Note that we neglect the propagation time of the light 
within the structure compared to particle motion. 

Using simple arguments based on the Maxwell stress tensor the momentary held then 
determines the optical force on the j-th particle along the x-axis [18, 21]: 

Fi ='j + \B,? - \Cj? - . ( 2 ) 

We will use this optical force then as ingredient to obtain Newton’s equations for the par¬ 
ticle motion, where we add damping describing additional light or background gas induced 
friction. In previous work we and others have seen that under quite general conditions 
stationary equilibrium positions can be found, where all particles are force free and feel a 
restoring force against small shifts. However, in contrast to a standard ID optical lattice, 
the absence of a force does not enforce an unaltered propagation of the light [18] but just 
guarantees a balance between internal and external light scattering [9]. Following Ref. [18] 
and including a linear friction term proportional to /r the equations of motion for the point 
particles of mass m then read: 
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mxj = —fiXj + Fj{xi,... ,X]si). (3) 

Stationary structures require ... ,xp^) = 0 for all j. Here these are not simple equidis¬ 

tant lattices, but rather complex particle arrangements with regular central parts including 
the formation of Bragg like reflectors at both ends of the equilibrium particle distribution 
denoted as x^p [9]. 


III. COLLECTIVE EXCITATIONS AROUND EQUILIBRIUM POSITIONS 

When a particle is weakly displaced from its equilibrium Xj{t) = it feels a 

restoring force and at the same time induces a global light field modihcation generating 
forces on all other particles. Once set free the perturbed particles start a correlated oscilla¬ 
tory motion around their equilibrium positions x^^\ In the following we study the spatial 
properties and eigenfrequencies of such small amplitude particle oscillations analogous to 
phonons of a lattice. Note that the existence of a stationary stable equilibrium conhgura- 
tion does not guarantee stability with respect to even very small such collective oscillations, 
which in general can grow in amplitude and disintegrate the particle array. 

A Taylor expansion of the forces gives the following set of coupled differential equations 

N 

m'ij = -fiij + Dji^i, (4) 

i=i 

with nonlocal coupling matrix 


Dji = -^Fj{xi = x^^\... ,xn = = lim ^Fj{x^ = x^°^ + v = 1,... ,N). (5) 

Note that these coupling constants would be an ideal basis for a coupled oscillator model 
in the ultracold gas limit. For the classical point particle model here we simply use the 
ansatz ^ in equations (4) and solve for the eigensystem of D, where the oscillation 

frequencies can be calculated from the eigenvalues via 

i/i ± a/— — 4mA,y 


oo„ = 


2m 


yielding the linearized solutions: 

F = (A jC 
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The eigenfrequencies include a real part describing damping or antidamping depending on 
its sign. If positive, it leads to an amplification of the oscillation amplitude and the particles 
can not form a stable stationary conhguration. Hence, to ensure stable configurations the 
eigenvalues have to fulhl the condition [18]: 

m(9(A,))2 < (8) 

Note that added damping can restore the stability only as long as all real parts of the 
eigenvalues are negative. Without external friction, /r = 0, Eq. (7) simplihes to 


= HjCW m -)- BjC W m ^9) 

and as long as we have no imaginary parts of the eigenvalues A^^ and the real parts are 
negative the particles are simply oscillating around their equilibrium position. 


A. Negligible mode coupling limit C = 0 


In the special case when particle mode coupling is very weak but we have a strong 
transverse pump held applied, one can set C = 0 while still keeping a nonzero scattering into 
the hber, thus r] ^ 0. Here the scattering of the held by the particles within the hber is 
neglected and any two particles interact equivalently. As predicted before [7] and also found 
by us in Ref. [9], in this limit the particles in a stationary state are equidistantly distributed 
and the transfer matrix can easily be calculated explicitly: 

Mtm = (Mbs ■ ■ Mbs = ■ Mbs- ( 10 ) 


Here the collectively scattered held intensities are symmetric to right and left and read 



and for the force on the j-th of N particles we get: 


F,= 


Pri COS {Nkd/2) sin ((2j — N — l)kd/2) 


( 12 ) 


sin {kd/2) 

with = Ij^/c the radiation pressure resulting from the transverse pump held. 

This force vanishes for diherent lattice constants d = {2n — 1)X/{2N) with n G N giving 


6 



potential stationary solutions. We will now calculate the linearized coupling matrices D for a 
small displacement of the /-th particle within such a solution. For the amplitude perturbation 
at the j-th cloud from displacing the /-th cloud hy ^ = ek we have to distinguish between 
the two cases if it is at the right or left side of the displaced /-th cloud: 
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with the perturbation b. For j < I the amplitudes can be calculated as follows 


(13) 
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Let us hrst look at the stability of the solution with the maximal possible lattice distance 
n = N, i.e. d = {2N — 1)X/{2N), where the force on the j-th particle induced by the 
perturbation of the /-th particle reads: 


^ , 2iV - 1, 

Fj { d = -A 

^ ' 2N 


and so Dji = ^Fj is: 


(15) 


Dji = P^k (^sin - 5ji cot ^ . (16) 

D is a symmetric matrix with all eigenvalues real and as it is a circulant matrix, they 
and the eigenvectors can easily be calculated to give: 


A.= 



/ 2 Tfi(iy-l)(j-l) 

(e N 


N 

J=1 ’ 


( 17 ) 


As we have seen in Eq. (9), stably oscillating configurations require real and negative 
eigenvalues. Hence this conhguration is stable and we hud N phonon modes with frequencies 
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(18) 


where we have dehned a; 2 ,o = y cis binding oscillation frequency for two particles at 
(^ = 0. As expected from translation invariance of the system, the lowest v = 1-mode 
corresponding to the center of mass motion has a zero eigenfrequency. Interestingly our result 
exactly coincides with the prediction based on nonlocal fiber enhanced resonant dipole-dipole 
interaction obtained before in Ref. [7]. In the limit of large particle numbers we hnd: 


ojy ^ 2(z/ - l)a;2,oi 


N 


(19) 


7r(3 -|- 4z/(z/ — 2)) 

Consequently the particle binding frequency grows with for large particle numbers 
stiffening the structure with particle number. Fig. 2 shows the phonon frequency as a 
function of mode number for different particle numbers. It demonstrates that the second 
and the A^-th mode posses higher frequencies than all the other modes which are almost 
degenerate. 

It is now interesting to see how the atoms interact with respect to their distance. For 
this we insert Eq. (5) into Eq. (4) and calculate the perturbation induced forces explicitly: 


( 5 ^ sin 

\¥j 
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exhibiting a coupling between the oscillators proportional to sin 


N 


( 20 ) 


depending on the 


distance \j — l\ and particle number N. Maximum coupling occurs between center and 
boundary particles at = 2|j — / | and it vanishes for \j — l\ N. 

Let us now insert other zero force solutions at smaller distances dn = into Eq. (5). 

Here we can explicitly calculate the eigenvalues in the limit of large particle numbers 


A - p k 

^’>''4(1/- 1 ) 2 -( 2 n - 1 ) 2 ' 


( 21 ) 


As according to Eq. (8) fully stable configurations are only possible if all eigenvalues are 
negative, we see that the hrst case, n = N, was the only stable conhguration with a distance 
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Figure 2: Phonon frequency as a function of the phonon mode i/ for C = 0 beginning with v = 2. 
The blue line corresponds to = 3, the red line to N = 5, the green line to N = 10, the violet 
line to iV = 20, the orange line to A" = 40 and the grey line to = 50. 

below A. An example how the particles evolve over time is plotted in Fig. 3, where we 
compare stable configurations setting n = N (Fig. (a)) with unstable ones setting n = N — 1 
(Fig. (b)). Note that starting from a stable configuration we obtain correlated oscillations 
of particles and helds around the stable position with a slow nonlocal excitation transfer 
between the sites. At negligible coupling, C = 0, one finds neither damping nor amplification 
of the oscillation modes but we still see a higher light intensity conhned at the center of 
the structure. Starting from an unstable zero-force conhguration the particles keep their 
position only for a short time until two particles collapse and the other particles reorder. 
Nevertheless, also this new conhguration is not stable and the whole order disintegrates. 

B. Collective dynamics for finite particle field interaction parameter C, 

Particles within the hber mode will not only couple pump light in and out of the hber, 
but also scatter it to the opposite propagation direction, into free space or absorb it. This 
is effectively described by a non-zero complex value of C [18]. As a hrst consequence this 
immediately leads to spatially inhomogeneous helds in the mode and thus non-equidistant 
stationary conhgurations, which, in general, do not allow for analytic treatment. Although 
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Position Xi[X] 


Figure 3: Ten particles evolving over time with only transverse pump and ^ = 0. In Fig. (a) 
we have chosen d = {n = N) as initial condition slightly perturbed by ^initial = 0.1 A, 

while in Fig. (b) the particles are initially positioned at an unstable zero force distance d = "^^^A 
{n = N — 1). r 2 ,o is defined as r 2 ,o = 

the equation for the force (2) still looks very simple, its explicit form gets already complicated 
for more than three particles and even the zero-force points can not be found analytically. 
Hence to still get some analytical insights we hrst study systems of only two and three 
particles. 


1. Two particles 


For two particles with no helds injected along the hber, the light scattered by one particle 
only interferes with the light scattered by the second particle. Translation and spatial 
inversion invariance of the system here implies a vanishing sum of the forces on the two 
particles as long as we neglect asymmetric chiral coupling [22]. The force as function of the 
real part Cr and imaginary part Q then depends only on particle distance d [9]: 


Fi — —F2 


Fj,|l — iCP cos{kd) 

4 (ICP + 0) cos2(f) + 2C.sin(A;d) + 1 ’ 


( 22 ) 
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and the zero force distances do = f ^ are independent on (. Interestingly 

even for strongly absorbing particles as e.g. gold nano-particles Q » (r, we thus obtain 
the same force free binding distances [4], Intuitively this can be understood from the 7r/2 
relative phase shift of the two scattered helds at this distance preventing interference. To 
study the stability and strength of binding we calculate the coupling matrix D; 


^ Pykl^ - K? (sin(kd) (2|CP + 2^ + 1) + 2(r) 

(2 (cos(/cd) -h 1) (ICP + (i) + 2sin(kd)(r + 1)^’ 


(23) 


and its eigenvalues. Inserting the zero force distances d = (^ -|- n) A the nonzero eigenvalue 
reads 

^ 2kP,\l-tC\^ 

^ l + 2(|CP + 0-C.)’ 

and for d = (| -|-n) A we get the same eigenvalue with opposite sign. As the non-zero 
eigenvalue for the latter case is real and negative, these are stable conhgurations. The 
phonon frequency characterizing the optical binding strength then reads: 



^ — ^2,0] 




1 + 2(|CP + Ci ~ Cr) 


(25) 


which for large values of ( stays hnite and converges to below the noninteracting value. 
Nevertheless for real ( U 2 reaches a maximum of 002 ^^/3 + at C — 0.618 (Fig. 5(a)), while 
we hnd a minimum coupling frequency for blue detuning at C — —1.618 of U 2 = 1 ^ 2,0 \/3 — 
(Fig. 5(b)). The dynamics in these two extreme cases is shown in hgure 5. Note that in 
both cases we can observe that the particles and helds oscillate and are periodically pushed 
apart whenever the intensity of the light trapped between them is maximal. Note that even 
in the blue detuned case for low held seeking particles we hnd stable trapping but much 
slower oscillations, i.e. weaker conhnement. 


2. Three particle dynamics 

Let us now add a third particle. Due to symmetry we will restrict ourselves to equidistant 
conhgurations and calculate the total transfer matrix. Again the amplitudes of the electric 
held to the left and right of the particles yield the forces acting on the particles, where as 
a consequence of symmetry the force on the middle particle is zero and the remaining two 
sum up to zero, i.e.: 
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Coupling constant (r 


Coni^ling constant Q 


Figure 4: Fig. (a) shows the phonon frequency oji^ as a function of the real part of the coupling 
constant Cr- The blue line corresponds to (i = 0, the red line to Ci = § and the green line to Ci = 5 - 
Fig. (b) shows the phonon frequency oji^ as a function of the imaginary part of the coupling constant 
(i. The blue line corresponds to Cr = 0, the red line to Cr = g and the green line to Cr = 5 - 
We can see that the function goes to for large values of C,. For the real part of C we can find a 
maximum for > 0 . 


Fi = —F 3 =~ Pr, {cos{kd) + cos{2kd) — ((2sm{2kd) + sm{3kd) + sin(4fc(i))) + 0[C]^, 

(26) 

F 2 = 0 

Here the general expressions are complex so that we only present some special cases to 
exhibit the stability of the particles. For real C = | only one stable equidistant configuration 
exists with di = d 2 — (0.8276 + ? 7 ,)A. Adding an imaginary part, ( = the stationary 
distances are enlarged to di = d 2 — (0.8305 + n)X. As shown in Fig. 6 the trajectories of 
the particles here show surprising features. 

Although the stationary state is symmetric with respect to the middle particle, a small 
perturbation from this equilibrium induces an intricate superposition of two oscillations, 
namely fast relative oscillations superimposed on a much slower and large amplitude a center 
of mass oscillation. As the system possesses translational invariance implying momentum 
conservation, it could be expected in principle to not allow any center of mass oscillations. 
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Figure 5: Two particles evolving over time with only transverse pump oscillating after a small 
perturbation with maximal frequency at C = 0.618 (a), and minimal frequency at C = —1.618 (b). 

This is at first sight also conhrmed by the oscillation eigenvalues {0, —3.83P^/c, —3.51P^/c} 
calculated for real C = 1/9 with eigenvectors 
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The eigenvalue Ai corresponding to the center of mass oscillation zi is indeed zero. However, 
we can see that the eigenvector Z 2 , corresponding to the case where two outer particles 
oscillate together against the middle particle includes a center of mass oscillation, as the 
middle particle moves much more then the sum of the outer two. This is also clearly 
visible in hgure 6 . At this point we have to recall that we have a strongly coupled atom 
held system where the pump laser constitutes an external energy and momentum reservoir, 
which can provide or accept momentum from the particles and conservation is only true for 
the combined atom held system. The asymmetric oscillation Z 2 thus periodically channels 
momentum in and out of the light held inducing these unexpected center of mass oscillations. 
The dynamics gets even more complicated if we allow for light absorption and investigate 
the eigensystem for an imaginary coupling constant, ( = We hnd the eigenvalues 

{0, —2.98P,jfc, —2.89Pj,A;} with eigenvectors 
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Figure 6 : Three particles evolving over time with only transverse pump, C = 5 (s-)) C = (b) 

and a small perturbation ^initial = 0.1 A. 
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(28) 


Again we observe a large amplitude center of mass oscillation due to the non inversion sym¬ 
metric initial perturbations which lead to even amplihed center of mass motion. Hence we 
see that while we have stationary equilibrium conhguration, the system tends to disintegrate 
due to dynamical instability when no external friction force is provided. 


IV. NUMERICAL STUDIES OF THE DYNAMICS OF LARGER ENSEMBLES 

While the momentary light held for a given conhguration of many particles can still be 
found by simple matrix multiplication, the explicit expression for the forces soon becomes 
prohibitively complicated and their common zeros cannot be analytically obtained any more. 
Numerically, however, this still can be done efficiently even for hundreds of particles and 
the solutions of the corresponding dynamical equations of motion are still possible. In order 
to capture typical features of the resulting many particle dynamics but still obtain readable 
plots, here we study the case of ten particles for various coupling strengths and damping. 
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Figure 7: Ten particles evolving over time with only transverse pump, C = g (a), C = ^ (b) 
and a small perturbation of (initial = O.OOIA. Fig. (c) shows the oscillation of the particles. Here 
we magnify the trajectories of Fig. (a) and plotted a range of O.OIA. In the second case the 
particles, hrst, start to oscillate around their equilibrium positions, but the perturbation leads to 
antidamping, which, in the end, destroys the configuration. The eigenvalues of the hrst case are all 
negative and real, while the eigenvalues of the second case acquire positive and negative imaginary 
parts. 


Fig. 7 compares the time evolution for real polarizability ( (Fig. 7(a)) with the case of partial 
absorption, where ( also possesses a small imaginary part (Fig. 7(b)). In the hrst case a 
small perturbation simply induces small correlated oscillations in the vicinity of local held 
maxima with some slow energy exchange between the particles. In the second case with 
absorption, the particles start to oscillate hrst, but these oscillations grow until one reaches 
the nonlinear regime leading to a complete destruction of the order. As shown in a recent 
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paper [9] one can still find a stable configuration for this case in the over damped limit. The 
mathematical origin for this instability is that due to the imaginary part of ( the eigenvalues 
also acquire imaginary parts. This instability still can be suppressed up to a limiting value 
by introducing an additional friction coefficient as in Eq. (8). 

We had demonstrated that in the weak-coupling-limit, C = 0) has stable oscillatory 
solutions without the need of damping. Introducing ( ^ 0 this behaviour is often lost 
and even for rather small perturbations the particle oscillations grow beyond the linearized 
regime. This leads to the collapse of the particle structure and they are expelled to the sides. 
It is now interesting to study the grade of instability in such a system, which turn out to be 
very parameter and size sensitive. To quantitatively show this we compute the eigenvalues 
of D and examine the real parts of them characterizing the magnitude of the corresponding 
instabilities. 

The maximum real part of all the nonzero eigenvalues of D is plotted in Fig. 8 as a 
function of real and imaginary part of ( for different system sizes N. If this value is negative, 
in principle one can always add enough extra damping to stabilize the conhguration, while 
a positive value implies dynamical instability in any case. The hgures show that for odd 
particle number we get a boundary line to the left of which at least one eigenvalue is positive 
and Eq. (8) thus predicts an instability, which can not be removed by damping. For negative 
real ( and even particle numbers (Fig. (b) and (d)) we hud stable conhgurations with 
negative eigenvalues but they are very very close to zero. In this case the system is very 
sensitive to small perturbations and needs strong extra damping for stabilization. Here 
in numerical simulations we often End unstable dynamics despite using a nominally stable 
parameter set, as validity of the linearized equations is very limited. 

An analogous phenomenon occurs for large particle numbers (Fig. (d)) and large values 
of Cr- In this case it is difficult to distinguish between stable conhgurations with eigenvalues 
close to zero and unstable configurations. We can observe that for large ( positive eigenvalues 
appear, especially for imaginary which prevent the formation of stable configurations. In 
general only for small and mostly real ( (blue region in Fig. 8 ) robust stable configurations 
can be expected. 
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Figure 8: Maximum real part of the nonzero eigenvalues as function of real and imaginary part of 
(. Fig. (a) corresponds to = 5 particles, (b) to = 10, (c) to A^ = 15 and (d) to A^ = 20. Black 
lines show the zero value contour lines. 

V. CONCLUSIONS 


Coherently illuminated particles in a ID trap order in strongly bound regular arrays 
induced by the interference of the scattered and pump light. If the scattered light is guided 
by nano-optical waveguide structures one hnds very strong binding and long range order 
with phonon frequencies growing approximately with the square root of the particle number. 
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In general the scatterers arrange in configurations, where the scattered light is resonantly 
confined within the structure to maximize the optical trapping potential for the particles. 
This can be seen as a controllable prototype ID implementation of an atom light crystal. 
In this regime absorption and scattering of the fields by the particles within the waveguide 
create a strongly nonlinear dynamical response of the system, which limits the size of stable 
arrays even in the presence of additional friction forces. This hints at the possibility to 
create self-ordered nano-cavity QED systems where the particles simultaneously form the 
resonator via atom mirrors [23] as well as the optical active system. As a next task towards 
a quantized field treatment one of course would need to estimate the magnitude of the self- 
ordered single photon coupling strength, which requires a sort of dynamic mode description. 
Generalizations to several frequencies, polarizations or higher order transverse modes should 
generate even more intriguing complex physics. In principle similar mechanisms should also 
be at work in planar 2D setups [13] leading to even stronger transverse light confinement 
and limiting the extension of optical matter. 
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